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^ : 1 Introduction 



We denote the set of trees with bounded maximum degree A by T^. Setting t„ = IT^"^], we 
introduce a generating function for these trees 

Let iiT be a given subtree. For a tree T^, we say that H occurs in T^ if there is a subtree 
in T^ isomorphic to H. Denote the number of occurrence of if in a tree T^ by t^. To count 
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the occurrences, we introduce a generating function in two variables as follows 



t{x,u)= Yl 



n t^ 
X U" . 



It can be simplified into 



t{x,u)= Y^ tn,kX'^v!', 



n>l,k>0 

where tn,k denotes the number of trees in 7^^ such that the number of occurrences of H is 
k. Note that t{x, 1) = t{x), i.e., t^ = Yl,k>o'^n,k- 

Furthermore, suppose that every tree in T^ is equally likely. We use X^iT^), or simply 
X„, to denote the number of occurrences of H in T^. Clearly, Xn is a random variable in 
T^. Then, the probability of Xn can be defined as 

Pr[X„ = k] = ^^. 

Without the maximum degree restriction, many results have been established for some 
special substructures. We further assume that, when H occurs in T^, the degrees of the 
internal vertices (vertices of degrees not equal to 1) should coincide with those of the corre- 
sponding vertices in the tree. In this case, H is called a pattern of the tree. Kok [S] showed 
that the number of occurrences of a pattern in the trees is with mean E{Xn) = {u + o(l))n 
and variance VariXn) = (a + o(l))n, and ?~ ^ " = is asymptotic to a distribution with 

yJVar{Xn) 

density {A + Bx^)e~^^ for some constants A,B,C > 0. Moreover, if the pattern is a star, 
then the number of occurrences of the pattern in a tree is exactly the number of vertices 
with a given degree with respect to the internal vertex of the star. It has been showed 

■y- rp/ ■v' ■! 

that for the number of vertices of a given degree Xn, ? " = is asymptotically normally 

distributed. We refer the readers to O [H] for more details. And, analogous results have 
been obtained for other classes of trees, such as simply generated trees, rooted trees, et al. 
(see |2], |S] and [5]). However, for the number of occurrences of a subtree H in trees without 
the maximum degree restriction, there is no such a similar result obtained. And, it seems to 
be much difficult. 

In this paper, we get that the number of occurrences of a subtree H in the planted trees 
and rooted trees with bounded maximum degree is also asymptotically normally distributed 
with mean and variance in 6(n). And for T^, we get a weak result which does not show that 
the distribution is also asymptotically normal. And, as an application, we use this result to 
estimate the Estrada index EE for these trees in T^, and give an explanation in theory to 
the approximate linear correlation between EE and the first Zagreb index [7] obtained by 
quantitative analysis. 

In this paper. Section 2 is devoted to a systematic treatment of the number of occurrences 
of a given subtree H. And in Section 3, we investigate the Estrada index for the trees in 

' n 



2 The number of occurrences of a given subtree 

In this section, we proceed to show that the number of occurrences of a subtree H in T^ is 
with mean {^h + o(l))n and variance [an + o{l))n for some constants ^h and an- And in 
the procedure of proof, we get the related results for planted trees and rooted trees. 

In what follows, we introduce some terminology and notation which will be used in the 
sequel. For the others not defined here, we refer to book [8]. 

Analogous to trees, we introduce generating functions for rooted trees and planted trees. 
Let TZ^ denote the set of rooted trees of order n with degrees bounded by an integer A. 
Setting r„ = |7?.^|, we have 

r{x) = y ^r„x" 

n>l 

and 

n>l,k>0 

where r„ ^ denotes the number of trees in TZ^ such that H occurs k times. A planted tree is 
formed by adding a vertex to the root of a rooted tree. The new vertex is called the plant, 
and we never count it in the sequel. Analogously, let V^ denote the set of planted trees of 
order n and bounded maximum degree A. Setting p„ = \V^\, we have 

p{x) = y^PnJ^" 

n>l 

and 



p{x, u) = ^ Pn,kX''u^, 
n>l,k>0 

where Pn,k denotes the number of trees in V^ such that H occurs k times. By the definition 
of planted trees, one can readily see that p{x, 1) = p{x) = r{x, 1) = r{x). 

Moreover, in [10], it has been showed that there exists a number xq such that 



p{x) = bi + h2\/xQ - X + h-i{xQ ~ x) ^ , (1) 

where 6i,62,&3 are some constants not equal to zero. Evidently, p{xq) = bi and for any 
\x\ < Xq, p{x) is convergent. For any A, xq < 1/2, particularly, if A = 4, xq ~ 0.3551817 
andp(xo) ^ 1.117421. 

Let ^'•^^^^(x) be the generating function of trees such that the degrees of the roots are 
not more than A — 1, while the degrees of the other vertices are still bounded by A. Then, 
we have (see [10]) 

p'-^-'Hxo) = I. (2) 

And, this fact will play an important role in the following proof. 



Let y(a;, u) = (yi(x, u), . . . , yN{x, u)Y be a column vector. We suppose that G{x, y, u) 
is an analytic function with non-negative Taylor coefficients. G{x, y, u) can be expanded as 

G{x, y,u) = ^ gn,kx''u''. 

n>l,k>0 

Let Xn denote a random variable with probability 

Pr[X„ = k] = ^, (3) 

9n 

where gn = J2k9n,k- First, we introduce a useful lemma OH]. 

Lemma 1. Let F{x,y,u) = {F^{x,y,u), . . . ,F'^{x,y,u))'^ be functions analytic around 
X = 0, y = {yi, . . . .yNY = 0, u = 0, with Taylor coefficients all are non-negative. Suppose 
F(0, y,u) = 0, F{x, 0,u) ^ 0, F^{x, y,u) ^ 0, and for some j, Fy_.y.{x, y,u) ^ 0. Further- 
more, assume that x = Xq together with y = y^ is a non-negative solution of the system of 
equations 

y = F{x, y, 1) (4) 

= det{I-Fy{x,y,l)) (5) 

inside the region of convergence of F, I is the unit matrix. Let y = {yi{x,u), . . . , yN{x, u)Y' 
denote the analytic solution of the system 

y = F{x, y, u) (6) 

with y{0,u) = 0. 

If the dependency graph Gf of the function system Equ. ^ is strongly connected, then 
there exist functions f{u) and gi{x,u), hi{x,u) (1 < i < N) which are analytic around 
X = Xq, u = 1, such that 



yi{x,u)=gi{x,u)-hi{x,u)Jl-—— (7) 

is analytically continued around u = 1, x = f{u) with arg{x — f{u)) ^ 0, where x = f{u) 
together with y = y{f{u),u) is the solution of the extended system 

y= F{x, y, u) (8) 

= det{I- Fy{x,y,u)). (9) 

Moreover, let G(x, y, u) he an analytic function with non-negative Taylor coefficients such 
that the point (xq, 2/(xo, 1), 1) is contained in the region of convergence. Finally, let X„ he the 
random variahle defined in Equ.^. Then the random variahle Xn is asymptotically normal 
with mean 

E{Xn) = fin + 0(1) (n -^ oo), 



and variance 

Var{Xn) = an + 0(1) {n -^ cx)) 

with jj, = t^. ' . 

Remark 1: We say that the dependency graph Gp of y = F(a;,y, m) is strongly connected 
if there is no subsystem of equations that can be solved independently from others. If Gp is 
strongly connected, then I — Fy(xo, Yo; 1) has rank A^ — 1. Suppose that v"^ is a vector with 
v^(I - Fy(xo, yo, 1)) = 0. Then, /x = ^ll^^(^Fi%'o,yo!i)) ■ ^^ ^^^^^ ^^^ readers to [21 i] for more 
details. 

Now, we concentrate on considering the generating function p{x,u). 

For the subtree H, we suppose the diameter of H is h. The depth of a vertex in a planted 
tree is the distance from the vertex to the root. The depth of a planted tree is the largest 
distance from the vertices to the root. We split up V^ into two sets Wq and W, which 
denote the trees with depth not more than h — 1 and the trees with depth greater than 
h — 1, respectively. We can see that if H occurs in the planted tree and the corresponding 
subtree takes the root, then the depth of the subtree is not more than h. Moreover, since we 
mainly consider the asymptotic number of subtree, the trees in Wq will contribute nothing 
to the coefficient of x"u'^ when n is large enough. Therefore, in this paper, we do not need 
to know the exactly expression of the generating function for the trees in Wo, and denote it 
by (f){x,u). Moreover, the counting function of a subset in Wo is also denoted by (f){x,u). In 
what follows, we shall see that this assumption is reasonable. Then, we focus on the trees 
in W. 

First, we introduce some conceptions. For a planted tree in W, the planted subtree 
formed by the vertices with depth not more than i is called i-depth subtree of this tree. 
Now, we split up W according to the h-depth subtree. That is, the trees in W having the 
same /i-depth subtree Wi form a subset Tii of W. Since the degrees of the vertices in W 
are bounded by A, there are finite number of different /i-depth subtrees, i.e., for some A^a, 
1 < i < A^A- Therefore, we obtain that 

p(x,m) = 0(x,m) + ^a^^_ft,(x,M), (10) 

where a^._/i(x, m) denotes the generating function of "Hj. 

To establish the functions system of a^-.h,(x,u), we need other functions a^'^/i_i(a;, u) as 
follows. For some {h — l)-depth subtree w'^, we denote the subset of the planted trees in W 
having w'^ by T-L'^. Note that wl ^ Ti'^. Then, the set of the planted trees having w^ consists of 
"H^ and w'^. If we use a^'.^h-ii^, u) to denote the generating function with respect to %[ Uw^, 
it follows that 

awr,h-i{x,u) = "^ ay,^^h{x,u)+w[{x,u). (11) 



There will appear an expression of the form Z{Sn, f{x,u)) (or /(x)), which is the sub- 
stitution of the counting series f{x,u) (or f{x)) into the cycle index Z{Sn) of the symmet- 
ric group Sn- This involves replacing each variable Si in Z{Sn) by f{x\u^) (or /(x*)). 
For instance, if n = 3, then Z{Ss) = {l/3\){sf + 3siS2 + 2S3), and Z^S^, f{x,u)) = 
(l/3!)(/(x, m)^ + 3/(x,m)/(x^,m^) + 2f{x^,u^)). We refer the readers to [8] for details. 

Note that a planted tree can be seen as a root attached by some planted subtrees. Em- 
ploying the classic Polya enumeration theorem, we have Z{Sj^i]p{x)) as the counting series 
of the planted trees whose roots have degree j, and the coefficient of a;^ in a; ■ Z{Sj^i]p{x)) 
is the number of planted trees with p vertices (see |8] p. 51-54). Therefore, 

A-l 

p{x) =x-^Z{Sj;p{x)), 

j=0 

and 

A-2 

p(^-^\x) = x-Y,Z{S,;p{x)). 

j=0 

By means of the same method, awi,h{x,u) can be expressed in a^'^ ^h-i{x , u) and </)(x, u). 
Suppose that the roots of the trees in "Hj have degree j, and each has j' planted subtrees 
with depth at least h — 1 attached to it. Clearly, / belongs to {1, ... ,j — 1}, and some of 
these subtrees may have the same w'^. Denote these different {h — l)-depth subtrees by {w'^} 
and suppose w'^ happens £s times. Evidently, ^is = j'- It follows that 



aw^^h[x,u) = X 



l[Z{i,-a^,^h-i) ■ (Pix,u) ■ u'^'-'^\ (1 < 2 < Na). (12) 



Here, (j){x, u) denotes the counting function of the other j — 1 —j' planted subtrees, since these 
subtrees belong to Wq. And the factor u''^^'"'f'^ serves to count the number of occurrences of 
H using the root of the new tree. In this case, all these vertices of the new tree corresponding 
the vertices of H have depth not more than h. And, since we know that the /i-depth subtree 
of the new tree is Wi, the number of occurrences taking the root can be calculated, that is, the 
upper index k{ls,(f)) can be calculated. Therefore, combining with Equ. Ollj) . the functions 
system of a«,i,/i(x, u) has been established. 

Now, we start to show that all the conditions of Lemma[T]hold for a^-^/i(x, u). For conve- 
nience, we still use F to denote the functions system. Set vector a(x, u) = {a^i^h, ■ ■ ■ , o,wn a)- 
We suppose aw.^h{x,u) = F*(x,a, u). Since p{x, 1) = p{x) and p{xo) = bi, one can see that 
awi,h{xo, 1) is convergent. So, xq and a(xo, 1) are inside the region of convergence of F. 
Apparently, the other conditions are easy to verify except for Equ.([5]). In what follows, we 
shall show that the sum Sa^_ ^ of every column of Fa(xo, a(xo, 1), 1) equals 1. Consequently, 
the equation det(I — Fa(xo, a(xo, 1), 1)) = holds. 

We consider the derivative on a^- ^h- If -F*(x,a, m) is not the function of a^- ^h, then 
F^ (x, a, u) will contribute nothing to the sum Sa^. ^- Thus, we just need to consider the 
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functions F*(a;,a, n) with some am/^^h-i having the term a^- ^h- In Equ. flT^ . if both a^'^ ^h-i 
and a^^ /i_i have the term a^i ,h, which imphes that the trees corresponding to a^/ ^h-i, 
ay^i^ ^h-i have the same {h — l)-depth subtree, then by the definition of ay^/^^h-i, we get that 
Oto^ ,h-i = ciw'^ ,h-i- Therefore, there exists exactly one product factor, say Z'(/s(,;au,^ ,h-i), 
in the expression of a^. ^h- 

Moreover, it is well-known that the partial derivative of Z{Sn', ■) enjoys (see [5]) 

d 

T; — Z[Sn', Sl, . . . , Sn) = Z{Sn-l', Si, ... , S„_i). (13) 

OS I 

For the planted tree, we have — g" (fn = Z{Sn^i;p{x,l)), which equals the generating 
function obtained by deleting one subtree from the root. Analogously, we have 



^. ._i)-0(x,m)-m'=(^-^). 



and it is exactly the new generating function produced by deleting one planted subtree of 
"H^y Uw^jj. Clearly, the root of the new planted tree is of degree j — 1. Particularly, if ig^ = I, 
after taking the derivative, the yielded function corresponds to the trees with roots of degree 
j — 1 such that every planted subtree does not belong to Ti'^^ U w'^^. Set u = 1. It follows 
that Sa^. ^{x, a(a;, 1), 1) is the generating function of the planted trees with roots of degree 
not more than A — 1, i.e., p^^~^\x, 1). Combining with the fact p^^~^\xo, 1) = 1, we obtain 
Sa^. hi^o, a(xo, 1), 1) = 1. Immediately, the Equ.([5]) 

det(I-Fa(xo,a(xo,l),l)) = 

follows. 

Employing Lemma [1], we have that a^.,ft,(x, u) is in the form of Equ.([7]), namely, for some 
f{u) and gwi,h{x,u), hyj.^h{x,u) which are analytic around x = Xq, u = 1, it follows that 



awi,h{x, u) = gwi,h{x, u) - /iu,,,/i(a;, u)Jl 



X 



fin) 



is analytically continued around u = 1, x = f{u) with arg(x — /(n)) ^ 0. From Equ. OlOl) . we 
can see that p{x, u) can be written into a function of a(x, u), and denote it by P(x, a(x, u),u). 
Clearly, all the coefficients of P(x, a(x, n), n) are non-negative. Therefore, p{x,u) is also in 
the form of Equ.([7]). Moreover, recalling Equ.([T]), we can see that /(I) = xq. Apply Lemma 
[T]to P{x,a{x,u),u), the following result is obtained. 

Theorem 2. For any given subtree H, the number Xn of occurrences of H in V^ is asymp- 
totical to be normal with mean E{Xn) = finn + 0(1) and variance Var{Xn) = cr^n + 0(1) 
for some constants fin o-nd a^. 

A rooted tree in TZ^ can also be seen as a root attached by some planted trees. That is, 
by the classic Polya enumeration theorem, analogous to Equ. (TT2!) . the generating function 



of TZ'^ is also a function in a(x, u). We denote the function by R{x, a(x, u),u), and r(x, u) = 
R{x,a{x,u),u). By means of the above analysis, it is not difficult to see that the Taylor 
coefficients of R{x,a{x,u),u) are non-negative. Thus, r{x,u) also has the form of Equ.©. 
And, apply Lemma [1] to R{x,a{x,u),u), the following result is obtained. 

Theorem 3. For any given subtree H, the number Xn of occurrences of H in TZ^ is asymp- 
totically normally distributed with mean E{Xn) = fJ^nn + 0{1) and variance Var{Xn) = 
o"]^n + 0(1) for some constants fin o,nd a^^. 

Remark 2: Since r(x,M) and p{x,u) correspond to the same function f{u), by Lemma [1] 
we can see that the means of Xn with respect to TZ^ and V^ are with the same constant 
fiH- Moreover, it has been showed that the sum of each column of Fa(xo, a(xo, 1), 1) equals 
1, then we have v-^ = (1, . . . , 1) such that v-^(I — Fy(xo, yo, 1)) = 0. Therefore, it is easy to 
see that fin is positive by Remark 1. 

In what follows, we investigate the generating function of trees. Two edges in a tree are 
similar, if they are the same under some automorphism of the tree. To join two planted trees 
is to connect the two roots with a new edge and get rid of the two plants. If the two panted 
trees are the same, we say that the new edge is symmetric. Then, we have the following 
lemma due to flOl. 



Lemma 4. For any tree, the number of rooted trees corresponding to this tree minus the 
number of nonsimilar edges (except for the symmetric edge) is the number 1. 

Note that, if we delete any one edge from a similar set in a tree, the yielded trees are 
the same two trees. Hence, different pairs of planted trees correspond to nonsimilar edges. 
Now, we have 

t(x, u) =r(x, m) - - ( ^ a^^^^h{x, u)a^^^^h{x, u) ■ u''^'"'^''"^^'^ 

l<ii,i2<NA 

+ 1 Yl a^.A^\u').u'^---'\ (14) 

l<i<NA 

where k{wij^,Wi^) serves to count the subtrees taking vertices both in Wi-^ and Wi^. Conse- 
quently, we obtain that t{x,u) is also in the form of Equ.([7]), i.e., there exist some functions 
g{x,u), h{x,u) which are analytic around x = Xq, u = 1, such that 



t{x,u) = g{x,u) — /i(x,m). /l 



is analytically continued around u = 1, x = f{u) with arg(x — f{u)) ^ 0. Here, we could 
not show that t(x, u) likes P(x, a(x, u),u) and R{x, a{x, u),u) that have non-negative Taylor 
coefficients, so Lemma [1] fails in this case. However, we can use the following result due to 
[9] to get a weak result for t(x, u). 
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Lemma 5. Suppose that t{x, u) has the form 



t{x,u) = g{x,u) — h{x,u),ll 



where ^{x,u), h{x,u) and f{u) are analytic functions around x = /(I) andu = 1 that satisfy 
h{f{l), 1) = 0, hxifil), 1) ^ 0, /(I) > and /'(I) < 0. Furthermore, x = f{u) is the only 
singularity on the cycle \x\ = \f{u)\ foru is close to 1. Suppose that X„ is defined as Equ.^ 
to y{x, u) . Then, E{Xn) = (/i + o{l))n and Var{Xn) = (a + o{l))n, where /i = — /'(1)//(1) 
anda = ^'' + ^Ji-nl)/f{l). 

Remark 3: This result does not tell us that the limiting distribution is asymptotically 
normal. If h{f{l), 1) ^ 0, this lemma is trivial by Lemma [H and if h{f{u),u) = 0, we can 
still get that the limiting distribution is normal by further analysis (see [5]). 

For t{x), it has been obtained that [TU] 

t{x) = Co + Ci{xo - x) + C2{xo - x)^/^ H , 

where Cq, Ci, C2 are some constants not equal to 0. Combining with the fact t{x, 1) = t{x), we 
can see that /i(/(l), 1) = and /i^.(/(l), 1) ^ 0. Moreover, the other conditions in Lemma [5] 
are easy to verify. Then, we formulate the following theorem. 

Theorem 6. Let Xn be the number of occurrences of a given subtree H in the trees ofT^. 
Then it follows that 

E{Xn) = ifiH + o{l))n 

and 

\/or(X„) = K + o(l))n, 

where fin o-nd a\j are some constants with respect to the subtree H . 

Following book [T], we will say that almost every (a.e.) graph in a random graph space 
Qn has a certain property Q if the probability Pr(Q) in Qn converges to 1 as n tends to 
infinity. Occasionally, we shall write almost all instead of almost every. 

By Chebyshev inequality one can get that 



?,[\X^-E(X^)\>n^"]<'^^^^0 



„,„ , V. as n — > oo. 



Therefore, for any subtree H^ X„ = [jih + o{\))n a.e. in T^. Then, an immediate conse- 
quence is the following. 

Corollary 7. For almost all trees in T^ , the number of occurrences of H equals {fiH+o{l))n. 



3 The Estrada index 

In this section, we explore the Estrada index for trees in T^. Let G be a simple graph with 
n vertices. The eigenvalues of the adjacency matrix of G are said to be the eigenvalues of 
G and to form the spectrum. Suppose that the eigenvalues of G are Aj, 1 < i < n. The 
Estrada index is defined as 

n 

EE{G) = J2^^^- 

4 = 1 

This index is invented in year 2000, and nowadays widely accepted and used in the information- 
theoretical and network-theoretical applications. And for this graph invariant, many results 
have been established. We refer the readers to a survey [3] for more details. 

Furthermore, for trees with n vertices, it has been showed that the path has the minimum 
Estrada index and the star has the maximum. And by quantitative analysis, there is an 
approximate linear correlation between EE and the first Zagreb index, i.e., "^df for trees. 
Denote "^df by D. That is, 

EE^aD + b, (15) 

where a and b are some constants. We refer the readers to [5] and [7]. 

In what follows, we shall get the estimate of EE for almost all trees in T^ and give an 
explanation to the correlation (IT5l) in theory. 

Denoting by Mk = Mk{G) = XlILi -^f ^^^ k-th spectral moment of G, and bearing in 
mind the power-series expansion of e^, we have 



EE{G) = J2 



Mk{G) 



k\ 

fc=0 

Note that Mi^{G) is equal to the number of closed walks of length k. For trees, one can 
readily see that 

fc=0 ^ '' 
Then, in a tree, the closed walk of length 2k forms a subtree with at most fc + 1 vertices. 
We have got that, for any given subtree, the number of occurrences of the subtree in T^ 
equals {^h + o{l))n a.e. Since there are finite different subtrees with at most k + 1 vertices, 
and each subtree corresponds to finite numbers of 2k closed walks, we can obtain that there 
exists a constant ^2k such that the number of 2k closed walk is {^2k + o{l))n a.e., namely, 

M2k = (/i2fc + o(l))n a.e. 
in T^. Moreover, we introduce a lemma due to Fiol and Garriga [6]. 
Lemma 8. For any graph G, M2k < 'Yll=i '^? 
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Recall that the degrees of a tree in T^ are bounded by A. So, X]r=i ^f'^ — '^'^''^ ^^'^ thus 

(2fc)! 



EE{T^) < e^n. Moreover, since ^^^g -^^ is convergent, for any positive number e, there 



exists an integer jo such that for any j > jo, J2k=j+i 72f§' ^ ^"^ Evidently, it is uniform for 
all the trees in T^. Therefore, we have 

fc=0 ^ ^ k=0 ^ ' 

Hence, we just contribute to consider the closed walks of length at most jo- 

For any integer j, we have ^^=0 tH^T — ^^' Therefore, X]fc=o 'mo!' ^^ convergent, and 
denote the limit by //a- It follows that 

(//A -e)n<^ — - = ^ -— < (/iA + o{\))n a.e. 

fc=0 ^ '' k=0 ^ '' 

Then, we have that (/iA — s)n < EE{T^) < (/ia + £)n a.e. Now, we can formulate the 
following theorem. 

Theorem 9. For any e > 0, the Estrada index of a tree in T^ enjoys 

{jj,A — £)n < EE{T^) < (/iA + e)n a.e., 
where fiA is some constant. 

If we suppose that the given subtree if is a path L of length 2, then there exists some 
constant ul such that in T^, the number of occurrences X„ of L is (ul + o(l))n a.e. In 
this case, it is easy to see that for each tree T^, Xn{T^) = ^^ (2) = \D{T^) — ?i + 1. 
Therefore, the value of Z) also enjoys {uD+o{l))n a.e. for some constant ud- Then, combining 
with Theorem [HI we can see that, for trees in T^, the correlation between EE and D is 
approximate to be linear. 
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